Non-Newtonian effects in the peristaltic flow of a Maxwell fluid 
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We analyzed the effect of viscoelasticity on the dynamics of fluids in porous media by studying 
the flow of a Maxwell fluid in a circular tube, in which the flow is induced by a wave traveling on 
the tube wall. The present study investigates novelties brought about into the classic peristaltic 
mechanism by inclusion of non- Newtonian effects that are important, for example, for hydrocarbons. 
This problem has numerous applications in various branches of science, including stimulation of fluid 
flow in porous media under the effect of elastic waves. We have found that in the extreme non- 
Newtonian regime there is a possibility of a fluid flow in the direction opposite to the propagation 
of the wave traveling on the tube wall. 
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>Y I. INTRODUCTION 

V 

^J | Investigation of flow dynamics of a fluid in a tube having circular cross-section, induced by a wave traveling on its 

FT* . • wall (boundary), has many applications in various branches of science. The physical mechanism of the flow induced 

by the traveling wave can be well understood and is known as the so-called peristaltic transport mechanism. This 

q , mechanism is a natural cause of motion of fluids in the body of living creatures, and it frequently occurs in the organs 

such as ureters, intestines and arterioles. Peristaltic pumping is also used in medical instruments such as heart-lung 

£>^j machine etc. [1]. 

Laboratory experiments have shown that an external sonic radiation can considerably increase the flow rate of a 

Oh liquid through a porous medium (Refs.[l,2] and references therein). Initially, the idea of flow stimulation via waves 

traveling on the flow boundary, in the context of porous media, has been proposed by Ganiev and collaborators [3] . 

,_^ ' They proposed that sonic radiation generates traveling waves on the pore walls in a porous medium. These waves, 

£> in turn, generate net flow of fluid via the peristaltic mechanism. Later, this problem has been studied in a number 

^sO ' of publications, where authors used different simplifying assumptions in order to solve the problem (see e.g. Ref.[4]). 

The most recent and general study of stimulation of fluid flow in porous media via peristaltic mechanism is presented 

in Ref.[l], which we will use as a starting point in order to include non-Newtonian effects into the peristaltic model. 

It is clear that a usual peristaltic mechanism discussed, e.g., in Ref.[l] can be used to describe the behavior of a 

classic Newtonian fluid; however, for example, oil and other hydrocarbons exhibit significant non-Newtonian behavior 

[5] . The aim of this paper is therefore to incorporate non-Newtonian effects into the classical peristaltic mechanism 

[1]. Thus, the present work formulates a realistic model of the peristaltic mechanism which is applicable to the 

non-Newtonian fluids (e.g. hydrocarbons) and not only to the Newtonian ones (e.g. ordinary water) which have been 

• T* , extensively investigated in the past [1] . 

It should be noted that there were similar studies in the past ([6] and references therein). However, the previous 

contributions discussed peristaltic mechanism for rheological equations other than the Maxwellian one. Thus, the 

Oh present study fills this gap in the literature. In addition, this study is motivated by the recent results of del Rio, 

de Haro and Whitaker [7] and Tsiklauri and Beresnev [8], who found novel effects, including the enhancement of a 

Maxwellian fluid flow in a tube that is subjected to an oscillatory pressure gradient. 



II. THE MODEL 

We consider an axisymmetric cylindrical tube (pore) of radius R and length L. We assume that elastic wave induces 
a traveling wave on the wall (boundary) of the tube with the displacement of the following form: 

W(z,t) = R + acos(—(z-ct)), (1) 

A 

where a is the amplitude of the traveling wave, while A and c are its wave-length and velocity, respectively. We note 
that z-axis of the (r,<f),z) cylindrical coordinate system is directed along the axis of the tube. 



The equations which govern the flow are the balance of mass 

dp 



and the momentum equation 



m , V • (pv) = 0, (2) 



"1" ' 

p-^+p(tV)v = -Vp-Vf, (3) 



where p, p and v are the fluid density, pressure and velocity, respectively; f represents the viscous stress tensor. We 
describe the viscoelastic properties of the fluid using the Maxwell's model [7], which assumes that 

im-^ = -MVv-|V-iT-r, (4) 

where /i is the viscosity coefficient and t m is the relaxation time. 
We further assume that the following equation of state holds 

-£=«, (5) 

where k is the compressibility of the fluid. We also assume that the fluid's velocity has only r and z components. 
We make use of "no-slip" boundary condition at the boundary of the tube, i.e. 

dW 
v r {W,z,t) = -z-, v z (W,z,t)=0. (6) 

at 

Eq.(4) can be re- written in the following form: 

Pi \ 

l + t m —)f = -pVv-^\7-v. (7) 

Further, we apply the operator (1 + t m d / dt) to the momentum equation (3) and eliminate f in it using Eq.(7): 

(8) 



d \ u ( d 



P-q7+P{VV)V 



The equations are made dimcnsionless by scaling the length by R and time by R/c. Also, we have introduced the 
following dimcnsionless variables (and have omitted the tilde sign in the latter equations): W — W/R, p — p/po, 
v r = v r /c, v z = v z /c, p = p/(pqc 2 ). Here, p is the regular (constant) density at the reference pressure p$. We have 
also introduced e = a/ R, a = 2irR/\, Re = pqcR/ p, \ = npoc 2 . 

Following Rcf.[l], we seek the solution of the governing equations in a form: 

p = Po + epi(r,z,t) + e 2 p 2 (r,z,t) + ..., 
v r = eui(r,z,t) + e 2 U2{r 1 z 1 t) + ..., 
v z = ev\(r,z,t) + e 2 v 2 (r,z,t) + ..., 

p = 1 + e/9i(r, z, t) + e 2 p 2 (r, z,t) + .... 

Then, doing a usual perturbative analysis using the latter expansions, we can obtain a closed set of governing 
equations for the first (e) and second (e 2 ) order. 

Further, following the authors of [1,8], we seek the solution of the liner problem in the form: 

Ul (r,z,t) = C/i(r)e M ( z -*) +C7 1 (r)e- M ( z - t \ 
Vl (r,z,t) ^V^ry^*-^ +V l (r)e- ia{z - t \ 



Pi 



(r,z,t) =Pi(r)e M ( z -*) + P^e^ 01 ^^, 



p 1 (r,z,t)^ X P 1 (r)e^ z -^+xPi(r)e- la(z - t \ 

Here and in the following equations the bar denotes a complex conjugate. 
On the other hand, we seek the second (e 2 ) order solution in the form: 

u 2 (r, z, t) = U 20 (r) + U 2 (r)e l2a(z ^ + U^e^ 2 "^, 
v 2 (r, z,t) = V 2Q (r) + V 2 {r)e aa(z ^ + V^e"^^, 



1>2 



(r,z,t) = P 20 (r) + P 2 (r)e i2a ^ z -^ +P 2 (r)e- l2a < z - t ), 



P'i 



(r,z,t) = D 2Q (r)+D 2 (r)e l2a(z ^+D 2 (r)e- l2a ^- t K 



The latter choice of solution is motivated by the fact that the peristaltic flow is essentially a non-linear (second 
order) effect [1], and adding a non-oscillatory term in the first order gives only trivial solution. Thus, we can add 
non-oscillatory terms, such as U 2 o(r), V 2 o(r), P 2 ${r), D 2 o(r), which do not cancel out in the solution after the time 
averaging over the period, only in the second and higher orders. 

In the first order by e we obtain: 

-(1 - iat m )P[ + -i- ( U'{ + ^ - % - a 2 U^\ +-L-^-(u{ + — + iaVi) = -ia(l - iat m )U u (9) 

He \ " r r z / ARe dr \ r J 

-»a(l - iat m )P[ + -±- \V(' + y - o?VA + ^- \U[ + ^ + iaVA = -*a(l - iat m )Vi, (10) 

U[ + — + iavA = iaxPi- (11) 

Here, the prime denotes a derivative with respect to r. 

Further, we re-write the system (9)-(ll) in the following form: 

-JP{ + (u[' + ^ - % - ^U.) = 0, (12) 

\ r r A 



I 



V{ 



-iPi --[V{' + -±- FVx = 0, (13) 

a \ r J 

where 

7 = (1 — iod m )Re — iax/3, (3 = a — ia(l — iat m )Re. (14) 

Note, that Eqs.(12)-(13) are similar to Eq.(3.11) from Ref.[l], except that 7 and f3 are modified by substitution 
Re — > (1 — icd m )Re. 

Repeating the analysis similar to the one from Rcf. [1] , we obtain the master equation for U\ (r) and find its general 
solution as 

U 1 (r) = C 1 I 1 (vr) + C 2 h(J3r), (15) 

where I\ is the modified Bessel function of the first kind of order 1; and C\ and C 2 are complex constants defined by 
_ a(3vil ({3) _ -a 3 il (v) 

^1 = r.r it r \t / a\ a — t t n\ T t M ; ^2 — 



2[a 2 I (v)h(/3) - pvl (/3)h( v )} ' 2[a*I o (v)I 1 (0) - MoWM] ' 



where 



2 2 

v = a 



(1 - x)(l - iat m )Re - (4/3)zax 
(1 — iat m )Re — (4/3)iax 



Here, Iq is the modified Bcssel function of the first kind of order 0. 
We also obtain the general solution for V\ (r) : 

iaC\ i[3C 2 , . 

Vi(r) = h{vr) H Io(pr). 

v a. 

The second-order solution V 2 o(t) can also be found in a way similar to the one used in Ref.[l]: 



^o(r) = D 2 - (1 - iat m )Re / [Vi(C)tMC) + Vi(C)tfi(C)]dC, 

Jr 

where D 2 is a constant defined by 

iaCi i/3 2 C 2 , . -iaCi i(/3 2 ~C 2 ) ,5, 

•°2 = 2~ I ^ V ) 2a~ ^' — 2 — ^ 2a J i(<?). 

The net dimensionless fluid flow rate Q can be calculated as [1]: 



Q(z,t) = 2ir 



e I vi(r,z,t)rdr + e 2 / v 2 (r 7 z, t)rdr + 0(e 3 ) 
o ^o 



In order to obtain the net flow rate averaged over one period of time, we have to calculate 

r 2n/a 



This time averaging yields 



a 

<Q>= — / Q(z,t)dt. 
27T Jo 



< Q >=2ire 2 / V 20 (r)rdr 
Jo 



or finally substituting the explicit form of V 2 q (r) we obtain for the dimensionless net flow rate 



< Q >= ire 2 



l 
D 2 - (1 - iat m )Re I r^V^U^r) + Ki(r)!7i(r)]dr 



n 



III. NUMERICAL RESULTS 



(16) 



In the previous section, we have shown that the inclusion of non-Newtonian effects into the classical peristaltic 
mechanism by using the Maxwell fluid model yields the following change: Re — > (1 — iat m )Re in all of the solutions. 

It is known that the viscoelastic fluids described by the Maxwell model have different flow regimes depending on 
the value of the parameter De = t v /t m , which is called the Deborah number [7]. In effect, Deborah number is a ratio 
of the characteristic time of viscous effects t v = pR 2 / ' [i to the relaxation time t m . As noted in Ref. [7], the value 
of the parameter De (which the authors of Ref. [7] actually call a) determines in which regime the system resides. 
Beyond a certain critical value (De c = 11.64), the system is dissipative, and conventional viscous effects dominate. 
On the other hand, for small De (De < De c ) the system exhibits viscoelastic behavior. 

A numerical code has been written to calculate < Q > according to Eq.(I6). In order to check the validity of our 
code, we run it for the parameters similar to the ones used by other authors. For instance, for e = 0.15, Re = 100.00, 
a = 0.20, x = 0, t m = we obtain < Q >— 0.2708706458, which is equal (if we keep 4 digits after the decimal point) 
to the result of the authors of Ref.[l] who actually obtained < Q >= 0.2709. 

Further, we have made several runs of our code for different values of the parameter t m . We note again that t m 
enters the equations because we have included non-Newtonian effects into our model. Eq.(16) will be identical to the 
similar Eq.(4.1) from Rcf.[l] if we set t m = in all our equations. 



The results of our calculation are presented in Fig.l, where we investigate the dependence of < Q > on the 
compressibility parameter \ for the various values of t m . In order to compare our results with the ones from Ref.[l], 
we have plotted < Q > for the following set of parameters: e = 0.001, Re = 10000.00, a = 0.001, t m = (solid line). 
We note that the curve is identical to the corresponding curve in Fig. 2 from Ref. [1]. This obviously corroborates the 
validity of our numerical code. Further, to investigate the dependence of the flow rate < Q > on t m , we perform 
the calculation for a few values of t m . When t m = 1.0, we notice no noticeable change in the plot as both curves 
coincide within the plotting accuracy. For t m = 100.00 (dashed curve with crosses), we notice slight deviation from the 
Newtonian limiting case (solid line), which translates into shifting the maximum towards larger x's. For t m — 1000.00 
(dash-dotted curve with asterisks), we notice further deviation from the Newtonian flow, which also translates into 
shifting the maximum towards larger %'s. However, for t m — 10000.00 (dashed curve with empty squares), we note 
much more drastic changes, including the absence of a maximum and rapid growth of < Q > in the considered interval 
of variation of compressibility parameter \. The observed pattern conforms to our expectation, since large t m means 
small De (De < De c ) and the system exhibits strong viscoelastic behavior. Note that t m is dimensionless and scaled 
by R/c. 

After the above discussion it is relevant to define quantitatively the transition point where the flow starts to exhibit 
(non-Newtonian) viscoelastic effects. It is known [7] that De = t v /t m — (pR 2 )/(nt m ). Now, using definition of 
Re = pcRj \i we can define critical value of t m as 

f Re\ R 
*** = [De-) ^ (17) 

In all our figures we have used Re = 10000.0. If we put the latter value of Re and the critical value of the Deborah 
number 11.64 [7] into Eq.(17) we obtain t m c = 859.11 (measured in units of R/c). Therefore, the values of t m greater 
than t m c (for a given Re) correspond to sub-critical (De < De c ) Deborah numbers for which viscoelastic effects are 
pronounced. 

In Fig. 2 we investigate the behavior of the flow rate < Q > on the parameter a, which is the tube radius measured 
in wave-lengths. Again, to check for the consistency of our numerical results with the ones from Ref.[l] and also 
investigate novelties brought about by the introduction of non-Newtonian effects (appearance of non-zero relaxation 
time t m ) into the model, we first plot < Q > versus a for the following set of parameters: e = 0.001, Re = 10000.00, 
X = 0.6, t m = 0. If we compare the solid curve in our Fig. 2 with the dashed curve in Fig. 3 of Ref.[l], we will note 
no difference, which again corroborates the validity of our numerical code. We then set t m to various non-zero values 
and investigate the changes introduced by non-Newtonian effects. As in Fig.l, we notice no change for t m = 1.0. For 
t m = 100.00 and t m = 1000.00, we notice that flow rate somewhat changes attaining lower values as a (radius of the 
tube) increases. 

We treat the latter, t m = 1000.00, case separately for the reason of an appearance of a novel effect of negative flow 
rates when the interval of variation of a is increased up to 0.05. Again, we expect that for large t m (t m > t m c — 
859.11), i.e. small De (De < De c ), the system should exhibit viscoelastic behavior. We note from Fig. 3 that for 
a > 0.035, < Q > becomes negative, i.e., we observe backflow. In fact, by doing parametric study we conclude that 
tmc is the critical value of t m above which we observe backflow. By increasing t m further, t m = 10000.00, we note from 
Fig. 4 that in this deeply non-Newtonian regime < Q > becomes highly oscillatory, but what is unusual again is that 
we observe the negative flow rates for certain values of a that is the tube radius measured in wave-lengths. Obviously, 
the negative < Q > means that flow occurs in the direction opposite to the direction of propagation of traveling wave 
on the tube wall. Oscillatory behavior (appearance of numerous of maxima in the behavior of a physical variable) in 
the deeply non-Newtonian regime is not new [8]. However, the flow of a fluid created by the peristaltic mechanism 
in the direction opposite to the direction of propagation of traveling wave is unusual and should be attributed to a 
complicated, non-linear form of the response of a Maxwell fluid to the stress exerted by the wave. 

IV. CONCLUSIONS 

In this paper, we investigated the dynamics of fluid flow in a tube with a circular cross-section, induced by a wave 
traveling on its wall (boundary). This problem has numerous applications in various branches of science, including 
stimulation of fluid flow in porous media under the effect of elastic waves. 

The present study investigates novelties brought about into the classic peristaltic mechanism by the inclusion of 
non-Newtonian effects based on the model of Maxwell fluid. 

We have found that in the extreme non-Newtonian regime there is a possibility of flow in the direction opposite to 
the propagation of the wave traveling on the tube wall. Somewhat similar effect is known as the acoustic streaming [9], 
in which an acoustic wave propagating in a tube induces a mean flow in the direction of propagation in the acoustic 



boundary layer, but in the opposite direction in the central part of the tube. The mean flow or acoustic streaming is 
caused by the presence of friction at the bounding surfaces of the tube. While fluid away from the neighborhood of 
a boundary vibrates irrotationally as the acoustic wave passes, fluid in close proximity to the boundary must vibrate 
rotationally to satisfy the no-slip condition on the tube wall. This deviation from inviscid, irrotational behavior 
provides an effective driving force known as the Reynolds stress. This effective force, because of it is quadratic rather 
than linear, has non-vanishing time-average tangential component to the tube wall that drives flow in the boundary 
layer. In the case considered in our paper instead of having acoustic wave propagating through the volume of the 
tube, we have a wave traveling of the tube wall, besides we have further complication of considering non-Newtonian 
(Maxwell) fluid (recall that the discovered effect is demonstrated for the case of non-Newtonian regime t m /t m c > 1.0). 
Similarly, the peristaltic flow itself is a second order non-linear effect. Therefore, the flow of a fluid created by the 
peristaltic mechanism in the direction opposite to the direction of propagation of traveling wave (i.e. back-flow) could 
be explained by a complicated, non-Newtonian, non-linear response of a Maxwell fluid to the stress exerted by the 
traveling wave. 
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Figure captions 

Fig. 1 Dimcnsionlcss flow rate < Q > as a function of compressibility parameter \. The parameters used are 
e = 0.001, Re = 10000.00, a = 0.001. t m = corresponds to the solid line, whereas t m = 100.00, 1000.00 and 
10000.00 correspond to the dashed curve with crosses, dash-dotted curve with asterisks, and dashed curve with empty 
squares, respectively. 

Fig. 2 Plot of dimcnsionless flow rate < Q > as a function of a. Here, e = 0.001, Re = 10000.00, \ = 0.6. t m = 
corresponds to the solid line, whereas t m = 100.00 and 1000.00 correspond to the dashed curve with crosses and 
dash-dotted curve with asterisks, respectively. 

Fig. 3 Plot of dimcnsionlcss flow rate < Q > as a function of a on a larger than in Fig. 2 interval of variation of a. 
Here, e = 0.001, Re = 10000.00, \ = 0-6, t m = 1000.00. 

Fig. 4 Plot of dimensionless flow rate < Q > as a function of a. Here, e = 0.001, Re = 10000.00, x = 0.6, 
t m = 10000.00. 



3.0 





2.5 


to 

1 

O 

H 

X 

A 


2.0 
1.5 


0> 
V 


1.0 




0.5 



1 


1 1 


1 ; 




/■*" ^x 






/X * 


\ s / *N 




/ 








\ v ' ~ ~ 


Jt , ' 


/' 




Z' 


.X 


v^ 


S, ' 


/ 


X v 




/' 




^- 


/"* 


,- H \* 


/ \ 


^' x 




<-*"" 


,,0 


\^ x . ~ 


-W" 


,Q''' 


T^ "- 




m ---H" 




i 


- -H" 

1 1 


1 



0.2 0.4 0.6 

X 



0.8 



I 

o 

H 

X 

A 

O* 
V 



1.6 - 

1.4 - 

1.2 - 

1.0 - 

0.8 - 

0.6 - 

0.4 - 

0.2 - 








\ h 




J I I L 



0.002 0.004 0.006 

a 



0.008 



0.01 



I 

o 

H 

X 

A 

0> 
V 




0.05 



a 



0.0008 
0.0006 - 
0.0004 - 
A 0.0002 - 

V o - 

-0.0002 - 
-0.0004 - 
-0.0006 




0.002 



0.004 



0.006 



0.008 



0.01 



Q 



